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Abstract 

The stability of strange quark matter is studied within the Nambu Jona- 
Lasinio model with three different parameter sets. The model Lagrangian 
contains 4-fermion (with and without vector interaction) and 6-fermion terms; 
the minimum energy per baryon number as a function of the strangeness 
fraction of the system is compared to the masses of hyperons having the same 
strangeness fraction, and coherently calculated in the same version of the 
model, and for the same parameter set. The results show that in none of the 
different parameter sets strangelets are stable, and in some cases a minimum 
in the energy per baryon does not even exist. 
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Strangelet detection in heavy ion collisions has been proposed long ago [1-5] as 
a signature of Quark Gluon Plasma (QGP) formation: it has been suggested that 
rather cold droplets of stable or metastable strange-quark matter may be distilled in 
heavy-ion collisions during the phase transition from a baryon-rich QGP to hadron 
matter: the prompt anti-kaon and pion emission from the surface of the fireball 
could rapidly cool the QGP, thus favouring the condensation into metastable or 
stable droplets of strange quark matter [4, 5] . 

The investigation of strangelet stability is therefore rather crucial: it is funda- 
mental to find out whether this system is more or less stable than hyperons, in order 
to understand which state is more likely to be produced in heavy ion collisions, ei- 
ther hyperons and strange mesons or strangelets. The properties and stability of 
strangelets have been discussed within the MIT bag model, a confined gas stabilized 
by the vacuum pressure B [4,6-9]; more recent calculations, employing different 
quark models [10-12], pointed out that strangelet stability is strongly model de- 
pendent. 

In this paper we want to discuss the properties of strange quark matter within the 
NJL model, which has already been successfully utilized in the past to study strange 
quark matter [13]; we consider homogeneous quark matter made up of u, d and s 
quarks; no /9-equilibrium is required: in the system there is a strangeness fraction 
Rs = Ps/p, P being the total baryon density of quarks and ps the baryonic density 
of strange quarks. Furthermore, electromagnetic interaction has been neglected, 
so that the minimum of the energy corresponds to an equal number of u and d 
quarks. We therefore consider the curves corresponding to the minimum energy 
per baryon number as a function of Rs and compare these curves to the hyperon 
masses evaluated in the same version of the NJL model and for the same parameter 
values, in order to understand whether strange matter is more stable than hyperons 
or viceversa. 

Many versions of the NJL model have been used in the past; in our calculation 
we use a three flavour Lagrangian of the form [14-17]: 
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In the above ■0=1 d is the quark field, ih = | is the mass 
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matrix, Ai . . . Ag are the Gell-Mann fiavour matrices, and Aq = -y^^- 

£(4) generates four-leg interaction vertices, while £(6) gives rise to six-leg, flavour- 
mixing, interaction vertices; G and Gy are parameters of the model with the dimen- 
sions of [L^] and is a parameter with the dimensions of [L^]. In the following we 
will consider both zero and nonzero values for Gy, in order to see the effects of a 
repulsive vector interaction on the stability of strangelets. In the limit = 0, Vi, 
the symmetries of the model Lagrangian are the following ones: 



Uvil) X SUv{3) X SUa{S) X SUc{3) 
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where, of course, SUc{3) is global; Ua{1) is broken by the existence of the axial 
anomaly. 

In the mean field approximation it is possible to evaluate the dynamical quark 
masses, generated by their interaction with the vacuum, which are given by the 
following gap equation: 



m* = mi- 2G{ipiipi) + 2K {ipjipj) {ipkipk) 
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where 
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is the chiral condensate for the i-flavour. Since the NJL model is not renormalizable, 
a three-dimensional regularization with an ultraviolet cut-off A is introduced. 
The energy density of the system in mean field approximation turns out to be: 
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In the above formula, the constant £0 is introduced in order to set the vacuum 
energy density equal to zero, and py. is the time component of the vector current: 
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py. = (ipi'jo'^i)- The dependence of the above formula on Rg and p can be easily 
found by recalling that: 

Ps = RsP (10) 

= Pv, = 3 p (11) 

Pv. = 3i?,p (12) 

in the above formulas, are the baryonic densities of quark i, which are related to 
PVi by: ^ 

Pi = ^PVi (13) 

and to the Fermi momenta by: 



= {Sn'psY^' (14) 
kF.,, - l^Pi^-Rs)] , (15) 



p being the total baryon number density in the system (p = N/V). In the above 
the colour degeneracy and baryon number 1/3 of the quarks have been taken into 
account^ . 

Prom the above formulas, the energy per baryon number turns out to be: 

Etot ^tot 

In order to investigate the influence of the parameter values on strangelet sta- 
bility, in this work three different parameter sets will be used: 



set 1 


set 2 


set 3 


m = rriu — ma — 5.5 MeV 
rUs = 140.7 MeV 
A = 602.3 MeV 
GA2 = 3.67 
KK^ = 12.36 


m = rriu — f^d — 5.5 MeV 
rUs = 135.7 MeV 
A = 631.4 MeV 
GA2 = 3.67 
KK^ = 9.29 


m = rriu — f^a — 3.6 MeV 
rris = 87 MeV 
A = 750 MeV 
GA^ = 3.67 
KK^ = 8.54 



These parameters have been employed in Refs. [17] (setl), [16] (set 2) and [18] (set 
3); in the first two cases, the current masses for the u and d quarks are fixed on the 
basis of isospin symmetry and of limits on the average rh = {rriu + ^d)/2 at 1 GeV 

^Colour degeneracy is taken into account, even if no explicit gluon field is present in the model 
Lagrangian. 
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scale, while the remaining four parameters are fitted to reproduce the masses of tt, 
K and rj' mesons, together with the pion-decay constant the third parameter 
set also reproduces fairly well these experimental values. 

We choose for Gy the two following values, in order to discuss the relevance of the 
repulsive vector interaction in strangelet formation: 

Gy = Gv^ O.SG (17) 

the second value being motivated, for example, in Ref. [14]. 

With the above parameter values, the effective quark masses in the vacuum and 
the chiral condensates turn out to be: 





set l(MeV) 


set 2(MeV) 


set 3(MeV) 




367.7 


335.5 


361 


771* 


549.5 


527.6 


501 




241.9 


246.7 


287 




257.7 


266.7 


306 



In Ref. [11] a comparison is made between the minimum energy per baryon with 
respect to p at fixed i?^, and the hyperon masses coherently calculated in the same 
model and for the same parameter values, in order to understand which state is more 
stable. In the same spirit of this work, we compare our curves to the theoretical 
hyperon masses evaluated in the NJL model and for the same parameter values: for 
the first two parameter sets and Gy = we used the hyperon masses evaluated in 
Ref. [19]; following the same techniques presented in Ref. [19], we also calculated 
the hyperon masses for parameter set 3 and for Gy = 0.5G: they are all shown in 
Table 1. In Fig. 1 we present our results for the three different parameter sets: the 
continuous lines correspond to Gy — 0.5G, the dashed lines to Gy = 0; concerning 
the hyperon masses, empty triangles correspond to Gy = and empty squares to 
Gy = 0.5G; full circles are instead the experimental hyperon masses. As it is evident 
from this figure, for all three parameter sets, with or without vector interactions, the 
curves corresponding to strangelets turn out to be well above the hyperon masses 
having the same strangeness fraction: independently of the parameter set used, 
and of the presence of the vector interaction, strangelets are, therefore, not stable 
in the NJL model. From the left panel we can see that with the parameter set 
one and a nonzero value for Gy, a minimum in the energy per baryon is present 
only up to Rs — 0.7, as already pointed out in Ref. [10]. In the central panel the 
results corresponding to parameter set 2 are shown, from which it is evident that 
the vector interaction further reduces the range of values of Rg corresponding to a 
minimum in the energy per baryon: in this case in fact a minimum is present only 
for 0.16 < -Rs < 0.56. Within parameter set 3 this range is even more reduced, as it 
appears from the third panel: in the case of Gy 7^ a minimum is present only for 
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0.21 < Rg < 0.45. With this parameter set, even for Gy = 0, the range of Rg which 
is compatible with the existence of a minimum is hmited to: R^ < 0.79. In order to 
get a feehng of the occurrence of these minima, we show the curves corresponding 
to the energy per baryon as a function of p for different values of Rg in Fig 2, for 
the three parameter sets and with Gy — 0.5G: full circles indicate local minima. 

Our analysis shows that the existence of stable or metastable strangelets is not 
supported by the NJL model: the curves corresponding to the minimum energy per 
baryon as a function of the strangeness fraction are always higher than the corre- 
sponding hyperon masses coherently calculated in the model for the same parameter 
values, and in some cases the minimum does not exist; these results seem therefore 
to indicate that hyperons are more likely to be produced in heavy ion collisions, since 
they are more stable. This fact confirms the model dependence of the strangelet 
stability, which could set serious challenges to the search for these objects in heavy 
ion collisions. 
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Baryon 


N 


A 


■^0 


Q- 


niexp (MeV) 


938.27 


1115.68 


1314.9 


1672.45 


niseti (MeV) 


970.86 


1096.34 


1274.51 


1493.10 


Gy^O m,et2 (MeV) 


928.03 


1067.12 


1261.68 


1486.26 


m,et3 (MeV) 


996.77 


1101.34 


1271.78 


1471.33 


ruseti (MeV) 


965.78 


1128.53 


1314.8 


1501.91 


Gy = Q.5G mset2 (MeV) 


918.78 


1096.21 


1287.73 


1479.47 


rusets (MeV) 


996.65 


1124.97 


1287.6 


1455.77 



Table 1: Experimental masses and theoretical masses of hyperons calculated in the 
NJL model. 
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Figure 2: Energy per baryon number as a function of p for different values of Rg'. the 
three panels correspond to the three different parameter sets and Gy — O.SG. Full 
circles indicate local minima. The dashed lines in the second and third panel indicate 
the first and last values of Rs corresponding to a minimum: 0.16 < Rs < 0.56, and 
0.21 <Rs< 0.45, respectively. 
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